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ABSTRACT 

We examined six exoplanet host stars with non-standard Hipparcos astrometric solution, which may be indicative of unrecognised 
orbital motion. Using Hipparcos intermediate astrometric data, we detected the astrometric orbit of HD 5388 at a significance level of 
99.4 % (2.7 cr). HD 5388 is a metal-deficient star and hosts a planet candidate with a minimum mass of 1.96 Mj discovered in 2010. We 
determined its orbit inclination to be / = 178.3^^^° and the corresponding mass of its companion HD 5388 b to be M2 = 69 + 20 Mj. 
The orbit is seen almost face-on and the companion mass lies at the upper end of the brown-dwarf mass range. A mass lower than 
13 Mj was excluded at the 3 cr-level. The astrometric motions of the five other stars had been investigated by other authors revealing 
two planetary companions, one stellar companion, and two statistically insignificant orbits. We conclude that HD 5388 b is not a planet 
but most likely a brown-dwarf companion. In addition, we find that the inclinations of the stellar rotation axis and the companion's 
orbital axis diff'er significantly. 
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1. Introduction 

For most extrasolar planets detected with the radial-velocity 
(RV) method, only their minimum masses are known because 
of their unknown orbit inclinations. Although the large num- 
ber of planets allows us to draw statistically sound conclusions 
about their mass distribution (Udry & Santos 2007; Howard et al. 
2010), it is highly desirable to derive this distribution without 
the inclination incertitude. Astrometric measurements can re- 
solve this ambiguity. For instance, Sahlmann et al. (2011) used 
Hipparcos data to analyse the astrometric motion of stars with 
potential brown-dwarf companions and found that about half of 
the candidates are low-mass stars. Here, we extend this work to 
stars with planet candidates and non-standard Hipparcos solu- 
tion, and discover the astrometric orbit of HD 5388, which hosts 
a recently announced planet candidate. 



2. Target selection and astrometric analysis 

In September 2010, the list of RV-planets at exoplanet . eu con- 
tained 461 entries around 389 stars. Only 286 host stars are in- 
cluded in the new Hipparcos reduction (F. van Leeuwen 2007). 
Many host stars of transiting exoplanets discovered in transit sur- 
veys (e.g., Kepler, CoRoT, HAT, WASP) and confirmed by RV 
are fainter than the completeness limit of Hipparcos. We selected 
the stars for which the new Hipparcos reduction found a non- 
standard solution and that were not flagged as member of a mul- 
tiple system, thus having solution types '1', '7', or '9'. Type '1' 
solutions are termed stochastic and adopted when the standard 
five-parameter solution (type '5') is not satisfactory and neither 
orbital nor acceleration models improve the solution in terms 
of The types '7' and '9' are given, when the model has to 
include proper-motion derivatives of first and second order to 
obtain a reasonable fit. Six stars satisfied these criteria and are 
listed in Table 4. 



HD5388 HARPS 




Fig. 1. Phase-folded radial velocities of HD 5388 as published 
by Santos et al. (2010). Red circles indicate the HARPS mea- 
surements and the solid line corresponds to the best-fit solution. 



The astrometric analysis was performed as described in 
Sahlmann et al. (2011), where a detailed description of the 
method can be found. We briefly recall the main elements of the 
analysis. Using the orbital parameters known from RV measure- 
ments, the intermediate astrometric data of the new Hipparcos 
reduction was fitted with a seven-parameter model depending 
on the inclination /, the longitude of the ascending node Q, the 
parallax m, and off'sets to the coordinates (Aa*, A^) and proper 
motions (A/^a* , A/^d) given in the published catalogue of F. van 
Leeuwen (2007). A two-dimensional search grid in / and Q de- 
fined starting values for a standard nonlinear -minimisation 
procedure identifying the global minimum. The uncertainties in 
the RV parameters are propagated to the astrometric solution by 
means of Monte Carlo simulations. The statistical significance 
of each astrometric orbit was determined by the distribution-free 



2 



J. Sahlmann et al.: HD 5388 b is a 69 Mjup companion instead of a planet 



Table 1. Orbital parameters of HD 5388 (Santos et al. 2010). 



Parameter 


Unit 


Value 


P 


(day) 


777.0 + 4.0 


e 




0.40 + 0.02 




(ms-^ 


41.7 + 1.6 


To 


(MJD) 


54570.0 + 9.0 


Ci) 


(deg) 


324.0 + 4.0 


M2 sin / 


(Mj) 


1.96 


(J sin / 


(mas) 


0.05 


rameters of the 


astrometric solution for P 


Parameter 


Unit 


Value 






191 


A^orb 




1.5 


O-A 


(mas) 


3.8 


/•-red 




1.02 


Significance 


(%) 


99.4 




(mas) 


-0.7 + 0.6 




(mas) 


0.9 + 0.6 


VT 


(mas) 


18.9 + 0.7 




(mas yr~ 


1) 0.5 + 0.6 




(mas yr~ 


1) -0.1+0.5 


Q 


(deg) 


298.0!^^.4 


^orbit 


(deg) 


178.3:JJ;^ 


a 


(mas) 


1.7 + 0.5 


M2 


(Mj) 


69.2 + 19.9 



permutation test employing 1000 pseudo orbits. Uncertainties 
in the solution parameters were derived by Monte Carlo simula- 
tions. This approach has proven to be reliable in detecting orbital 
signatures in the Hipparcos astrometric data and efficiently dis- 
tinguishing a significant orbit in the present low signal-to-noise 
ratio regime (Sahlmann et al. 201 1). 

3. The orbit of HD 5388 

HD 5388 (HIP 4311) is Hsted with a stochastic solution in the 
catalogue of the new Hipparcos reduction. Santos et al. (2010) 
used the HARPS spectrograph to discover a massive planet can- 
didate around this F6V-star. Figure 1 shows the radial velocities 
from the discovery paper. The companion has a minimum mass 
of M2sin/ = 1.96 Mj and orbits its host in an eccentric orbit 
(e = 0.4) with an orbital period ofP = 111 days. The time of 
periastron passage To, the longitude of periastron oj, the radial- 
velocity semi-amplitude Ki, and the minimum semimajor axis 
of the star's astrometric motion a sin / are given in Table 1. The 
residuals of the RV solution (3.3 ms"^) are larger than the av- 
erage measurement error (2.8 ms"^). By considering the star's 
spectral type and projected rotational velocity and analysing the 
bisector inverse slope, Santos et al. (2010) conclude that the ex- 
cess RV noise is insignificant. 

We detected the astrometric orbit with a significance of 
99.4 % (2.7 (t), which means that Hipparcos measured the stel- 
lar orbit induced by a companion with the orbital characteris- 
tics given in Table 1. The derived astrometric orbit with a = 
1.7 ± 0.5 milli-arcsec (mas) is small compared to the median 
single-measurement precision (Ta = 3.8 mas given by the new 
Hipparcos reduction. However, the large number of measure- 
ments A^mes = 191 covering more than one orbital revolution 
(NoY\y = 1.5) resulted in an eflTective signal-to-noise ratio of 
S/N = a- ((Ta/ V^mes)"^ = 6.2 allowing the detection to be pos- 
sible. With a reduced chi- square value of x^^^^ = 1-02 the final 
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Fig. 2. Top: Astrometric orbit of HD 5388 projected on the sky. 
North is up and east is left. The solid red line shows the model 
orbit, which is orientated clockwise, and open circles mark the 
individual Hipparcos measurements. Bottom: 0-C residuals for 
the normal points of the orbital solution (filled blue circles) and 
the standard five-parameter model without companion (crosses). 



seven-parameter fit is good and the derived oflTsets to positions, 
parallax, and proper motions are small. The values of Q and the 
inclination /orbit are well constrained at Q = 298^2?° ^orbit = 
178.3+^-4°. Using a stellar mass of 1.2 (Santos et al. 2010), 
we derived the companion mass to be M2 = 69.2 ± 19.9 Mj. The 
astrometric orbit of HD 5388 is shown in Fig. 2 and its char- 
acteristics are summarised in Table 2. Reflfert & Quirrenbach 
(2011) found an upper mass limit of 124 Mj for the companion 
of HD 5388, which is compatible with our result. 

Assuming that the orbit orientations are randomly oriented 
in space, the probability of measuring an inclination smaller than 
2° is only 0.06 %. However, the number of planets detected in 
RV-surveys reached several hundreds making our detection rea- 
sonably probable. HD 5388 is a metal-deficient star with [Fe/H] 
= -0.27 ± 0.02. The hydrogen-burning mass limit, commonly 
used to distinguish between stars and brown dwarfs, increases 
with decreasing metallicity and reaches 87 Mj (0.083 Mq) for 
[Fe/H] < -1 (Chabrier & Baraflfe 1997). Because of the moder- 
ate metal-deficiency of HD 5388, we used a limit of 80 Mj and 
found that the probability that the companion of HD 5388 has 



0.01 



0.008 



J. Sahlmann et al.: HD 5388 b is a 69 Mjup companion instead of a planet 3 

Table 3. Spin inclination, orbital inclination, and obliquity. 
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Fig. 3. Probability density functions of orbit inclination (/orbit, 
dark-grey histogram) and stellar spin axis inclination (/rot, light- 
grey histogram). The orbit inclination PDF is very sharply 
peaked with its maximum at 0.18 and is truncated for display 
clarity. The PDF ofi/^ = 

/orbit ~ /rot IS shown as dashcd line. The 

bin size is 0.2°. 



a mass below this limit is 71 %. The companion of HD 5388 is 
thus most likely a brown dwarf. This adds HD 5388 to the very 
few known Sun-like stars that have a brown-dwarf companion 
and a well-determined astrometric orbit. 



3.1. Spin-orbit alignment 

After determining the orbital inclination from astrometry, we 
compared it to the orientation of the stellar spin axis. The in- 
clination /rot of the stellar spin axis is defined with respect to 
the line of sight (/rot = 0° or 180° for a pole-on view) and 
can be derived from the spectroscopic estimate of sin /rot = 
4.2 km s"^ (Santos et al. 2010). The authors did not give an 
error bar for this measurement and we assumed a conservative 
uncertainty of 1 km s"^. On the basis of the activity indicator 
log 7?^^ determined from the HARPS spectra, we used the cal- 
ibration of Mamajek & Hillenbrand (2008) to derive the stel- 
lar rotation period P^ot- Assuming an eff'ective temperature of 
Teff = 6297 ± 32 (Santos et al. 2010), the apparent visual magni- 
tude my = 6.839 ± 0.001 (F. van Leeuwen 2007), and the paral- 
lax m given in Table 2, we derived the star's absolute magnitude, 
luminosity (L = 4.5 ± 0.3 L©), and radius (R = 1.8 + 0.1 7?©), us- 
ing standard formulae and Monte Carlo resampling. Bolometric 
corrections were computed using the Flower (1996) parame- 
ters given by Torres (2010). Using the stellar radius and rota- 
tion period, we derived the equatorial rotation velocity to be 
V = 10 + 3 kms"\ and then the inclination of the spin axis by 
calculating /rot = arcsin(L'sin/rot/^'). 

Because we cannot determine the star's sense of rotation, the 
angle /rot has a 180° ambiguity and we first considered the value 
that falls into the same quadrant as /orbit, i-^- a prograde configu- 
ration. The astrometric analysis yielded the distribution of /orbit, 
obtained from 100000 Monte Carlo simulations. To obtain the 
/rot-distribution, we performed 10^ Monte Carlo simulations. We 
finally compared these two distributions by drawing 2-10^ pairs 
of values [/orbit, /rot] and obtained the distribution of the orbit 
obliquity or spin-orbit angle i// by computing i// = /orbit - hot- 
Fi gure 3 shows the probability density functions (PDF) of /orbit, 
/rot, and i//. The /orbit-distribution is very narrow and peaks around 



Parameter 


Unit 


Value 


1 <T interval 


3 <T interval 






-5.00 + 0.02 






Prot 


(day) 


9.4 + 2.0 






hot 


(deg) 


154.8!?o^, 


(144.5, 162.9) 


(111.1,174.4) 


^orbit 


(deg) 


178.3:0-4 


(177.6, 178.7) 


(171.4,179.1) 




(deg) 


23.4!i0-3 


(15.2,33.7) 


(3.3,67.0) 



North 
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Fig. 4. Possible spin-orbit configurations. The astrometric mo- 
tion in the sky plane defined by the north and east axes is ob- 
served from location P. The orbit orientation vector (red arrow) 
is defined by the angles /orbit and Q. The stellar spin orientation 
vector can lie on either the dark-grey cone (prograde configu- 
rations) or the light-grey cone (retrograde configurations) with 
identical opening angles 2 • (180° - /rot) = 51°. 



178.3°, whereas the /rot-distribution is broad with its maximum 
at ~ 155° and the two distributions show a very small overlap. 
The median values and confidence intervals of /rot, /orbit, and 
the obliquity i// are given in Table 3. We found that i// is larger 
than 4.4° and 19.3° with a probability of 99.7% and 68.3 %, re- 
spectively, and its nominal value with 1 cr confidence intervals is 
i// = 23^g^° in prograde configuration, which indicates that the 
orientations of stellar spin and orbital axis diff'er substantially. In 
retrograde configuration, the obliquity would be i//' = 153^^^°. 

We note that these values of i// and i//' are lower limits, be- 
cause the ascending node Qrot of the spin axis is not constrained 
and we assumed that D = Qrot in their derivation. As illus- 
trated in Fig. 4, i// increases if the two ascending nodes Qrot 
and do not coincide, creating an additional uncertainty in i// 
of 2 • (180° - /orbit) = 3.4° in the prograde configuration. For 
a retrograde orbit, this uncertainty is larger at 51°. In summary. 



the obliquity is i// = 23^_l^^^^-^^" and i//' = 153^^^^^^^^° for pro- 
grade and retrograde orbits, respectively, where the additional 
uncertainties are denoted in square brackets. 

We performed these calculations using the Prot-calibration 
by Noyes et al. (1984), which yielded a value of 8.8 + 1.6 days 
and an obliquity ofi// = 22^^°, in agreement with the result ob- 
tained with the calibration of Mamajek & Hillenbrand (2008). 
The isochrones of Marigo et al. (2008)^ indicated that HD 5388 
is at the end of the core-H burning phase just past the turn-off 
point. Therefore, HD 5388 may be about to leave the main se- 
quence or even be slightly evolved. This would introduce ad- 
ditional uncertainties and possibly biases into our estimates of 



The interface is http://stev.oapd. inaf. it/cgi-bin/cmd. 
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Table 4. Stars with candidate RV-planets and non-standard 
Hipparcos astrometric solutions. 



Object 


HD 


HIP 


Sol. 


Sp.T. 








type 




HD 5388 


5388 


4311 


1 


F6V 


55 Cnc 


75732 


43587 


1 


G8V 


HD 81040 


81040 


46076 


7 


GOV 


HD 179949 


179949 


94645 


1 


F8V 


HD 195019 


195019 


100970 


7 


G3IV 


7 Cep 


222404 


116727 


1 


KlIV 



log 7?^^, P^ot, the stellar radius, and eventually /rot, and may 
therefore alter the above findings. 

4. The orbits of five additional selected stars 

In addition to HD 5388, we examined five other stars that satis- 
fied our selection criteria in Sect. 2. They are listed in Table 4 
and although their astrometric motions were previously studied 
in the literature, we describe our findings below: 

- 55 Cnc has five planetary companions announced by Fischer 
et al. (2008). The planet 55 Cnc d has a minimum mass of 
M2 sin / = 3.84 Mj. Using HST astrometry, Mc Arthur et al. 
(2004) derived an orbit inclination of 53 + 7° for the d-planet, 
an astrometric perturbation size of a = 1.9 ± 0.4 mas, and a 
mass of M2 = 4.9 ± 1.1 Mj. Compared to the measurement 
precision of (Ta = 1 .4 mas, the stellar orbit is large enough to 
have been measured by Hipparcos, which explains the non- 
standard solution type. Because the Hipparcos data only cov- 
ers 20 % of the 5200 day orbital period, we did not attempt 
to fit an astrometric orbit. 

- HD 81040 hosts a massive planet with M2sin/ = 6.86 Mj 
(Sozzetti et al. 2006). We did not detect orbital signature be- 
cause the obtained orbit has a very low significance of 68 %, 
in agreement with the result of Sozzetti et al. (2006). As in 
that work, we used the Hipparcos astrometric data to set an 
upper limit to the companion mass of 47.5 Mj, thus exclud- 
ing that the companion has a stellar nature. 

- HD 179949 hosts a planet candidate with M2 sin / = 0.84 Mj 
(Tinney et al. 2001). We did not detect an orbital signature 
because the obtained orbit has a very low significance of 
71 %, which confirms the result of Zucker & Mazeh (2001). 

- HD 195019 hosts a planet candidate with M2 sin / = 3.70 Mj 
(Wright et al. 2007). We detected the astrometric orbit with 
a high significance of 98 % and a very low inclination of 
0.3 + 0.1°, which renders a stellar companion and confirms 
the resuk of Zucker & Mazeh (2001). 

- y Cep (HD 222404) hosts a planet candidate with M2 sin / = 
1.7 Mj (Hatzes et al. 2003). We did not treat this system be- 
cause an exhaustive analysis, including the use of Hipparcos 
astrometry, was performed by Torres (2007), who set an up- 
per mass limit of 16.9 Mj to the planetary companion. 



5. Discussion 

Among the potential planets around the six stars selected in 
Sect. 2, one turned out to be a star (HD 195019 b), two are likely 
planets (55 Cnc d and y Cep b), and one is a brown-dwarf com- 
panion (HD 5388). The astrometric orbits of two stars were not 
detected with Hipparcos and the companions remain planetary 
candidates. The significant fraction of one third of non-planetary 



companions found in this small sample is a selection eff'ect and 
does not reflect the properties of the planetary population. That a 
large population of non-planetary companions is observed with 
low orbital inclinations, hence mistaken for planets, is statis- 
tically unlikely and would have been detected with Hipparcos 
(Zucker & Mazeh 2001). At the individual level, however, a 
planet detected by RV remains a planet candidate until its mass 
range is confined by complementary constraints derived for in- 
stance from a transit lightcurve, dynamical interactions, by direct 
imaging, or from astrometric data, as in the case of HD 5388. In 
the future, as the data of the GAIA astrometry satellite becomes 
available, we will be able to measure the astrometric orbits of 
many planetary systems, thus obtaining an accurate census of 
the planetary mass distribution. A few more planet candidates 
may then turn out to be brown-dwarf companions. 

We have detected the astrometric orbit of HD 5388, which 
harboured a planet candidate of minimum mass M2 sin / = 
1.96 Mj. We found an almost face-on orbit with an inclination of 
/orbit = 178.3^Q^°. Consequently, we determined the mass of the 
companion Hij 5388 b to be M2 = 69+20 Mj. A mass lower than 
13 Mj was excluded at the 3 cr confidence level. Thus, HD 5388 b 
can no longer be considered a planet. It is instead a close brown- 
dwarf companion, which is remarkable because at most 0.6 % of 
Sun-like stars have close brown-dwarf companions (Sahlmann 
et al. 2011). In particular, it is one of the very few brown-dwarf 
companions for which the astrometric orbit could be measured 
and whose orbital motion is therefore fully characterised, such 
as HD 38529 (Benedict et al. 2010). 

The measurement of the spin-orbit alignment in planetary 
and binary systems, e.g., Triaud et al. (2010) and Albrecht et al. 
(2011), gives insight into the formation of these systems. The 
most frequently applied technique uses spectroscopic measure- 
ments during the transit of eclipsing systems and only deter- 
mines the projected spin-orbit angle. Because it is restricted to 
eclipsing systems, it favours small- separation, short-period sys- 
tems (another method relies on interferometric imaging of the 
stellar surface, see Hinkley et al. 201 1). We determined the stel- 
lar spin axis inclination on the basis of the spectroscopic estimate 
of sin /rot and the rotation period calibration using log 7?^^. 
We found evidence of non-zero spin-orbit angles with numeri- 
cal values of <A = 23!g^^^^ '^^° and il/' = 153^^(5^^^^° for prograde 
and retrograde orbits, respectively. Thus, for both orientations of 
the unknown stellar spin, we found substantial non-coplanarity 
between the stellar and orbital axes. Spin-orbit angle measure- 
ments obtained from the determination of the astrometric orbit 
and the spectroscopic v sin / off'er the opportunity to access this 
information for non-eclipsing systems. 

Acknowledgements. J. S. thanks A. Triaud for his contributions to the discussion 
on spin-orbit angles. 
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